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We show the existance of the exact plaquette-ordered ground states of the Hubbard model in- 
cluding site-off-diagonal interactions in arbitrary dimensions, by decomposing the Hamiltonian 
as sum of products of projection operators for each spin sector. The obtatined exact ground 
states are interpreted as Neel ordered states on the dual lattices. We demonstrate this idea 
in the one-dimensional chain and higher-dimensional lattices such as the Kagome lattice, and 
determine parameter regions of the exact ground states. 
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1. Introduction 

The Hubbard model is one of the generic models to 
describe interacting electrons in narrow-band systems. 1 
This model has played central roles to study magnetism 
and superconductivity. In spite of its simplicity, it is dif- 
ficult to solve this model exactly except for one dimen- 
sion or some special cases. On the other hand, it has 
been considered many extensions of the Hubbard model. 
The on-site repulsion of this model is due to the ma- 
trix elements of the Coulomb interaction corresponding 
to the on-site Wannier states, and other matrix elements 
are neglected. Therefore, we consider effects of these ne- 
glected terms as site-off-diagonal interactions. 2 For these 
generalized models, exact results for ferromagnetism and 
superconducting states have been discussed. 3-9 

Few years ago the authors discussed a different type 
of exact ground state in a one-dimensional system that 
is "bond Neel" state, 10 using the decomposition of the 
Hamiltonian into projection operators. 11 Furthermore 
they suggested extensions of this argument to higher di- 
mensional systems. 10 Main aim of this paper is to demon- 
strate the existence of the exact ground states in the 
generalized Hubbard model in arbitrary dimensions, es- 
pecially that on the Kagome lattice where few results 
have been known for Hubbard- type models. 12 ' 13 

The rest of this paper is organized as follows: In sec. 2, 
we review the method to construct Hamiltonians with ex- 
act ground states. In sec. 3, we apply this method to the 
one-dimensional model discussed in rcf. 10, and obtain 
the extended version of the phase diagram. In sec. 4, we 
apply the analysis to the Kagome lattice. Finally, we give 
summary and discussion of the results. 

2. Method 

The strategy to construct a Hamiltonian with an exact 
ground state is the following way. 11 First, we consider a 
Hamiltonian that can be decomposed as sum of products 
of operators for different (spin) sectors as, 



where a denotes the position of one of the unit plaquc- 

t , 4 ) is an operator 
0. The expectation value of this 



ttes that cover the lattice. R^J (a 



satisfying [R^,R^] 
operator is nonnegative (Raa) > 0. This condition is re 
alized, if R^a} is given by a product of an operator and 
its Hermitian conjugate. Then the expectation value of 
the Hamiltonian is also nonnegative (H) > 0. 

Next, we introduce a trial wave function given by a 
direct product of up and down spin sectors, 



|*(4B)> = I*tCA)>®I*i(S)>. 



(2) 



where A and B denote two groups of plaquettes that 
cover the lattice satisfying AUB = {all lattice sites}. We 
require that the projection operators have the following 
conditions, 



rM\^(A))=R$\$ 1 (B)) = 0, 
where a £ A and (3 £ B. Note that 
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Therefore, the eigenvalue of the Hamiltonian for 
| ^(.4, B)) is always zero. Then, the lower bound and the 
upper bound of the energy coincide, so that \^>(A,B)) 
turns out to be one of the exact ground state of this 
system. 

The above argument is satisfied in corner sharing lat- 
tices with the bipartite structure in terms of their dual 
lattices. For example, in two dimensions, the Kagome 
(Checkerboard) lattice can be covered by two plaquettes 
for different sectors alternatively, as illustrated in fig. 1. 
These states can be regarded as the Neel ordering on 
the honeycomb (square) lattice. In three dimensions, the 
Pyrochlore lattice and the Garnet lattice 14 satisfy these 
conditions. 

3. One-dimensional chain 

First, we consider the one-dimensional generalized 
Hubbard model at half-filling and zero-magnetic field, 
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Fig. 1. Examples of lattice structure where generalized Hubbard 
models with exact plaquctte-ordered ground states can be con- 
structed: the one-dimensional chain, the Kagome, the Checker- 
board and the Pyrochlore lattices. The solid and the dashed pla- 
quettes denote those belong to the groups A and B, respectively. 

given by H = ~Y^ ia ft-i.i+i,cr with the local bond Hamilto- 
nian, 

+V\\n la n J<y + V±_n la n-j S 

■j-X Tij a (tIjct + n ja) H ^~ / t Tij a Tij a i , (5) 

a' 

where a is the opposite spin of tr, and periodic boundary 
conditions are assumed. The number of the nearest sites 
is z = 2. We have defined the hopping and the density 
operators as Tij a = cJ^c-^ + H.c, = c\ a c ia . Note that 
the bond-bond interaction (W) term can be rewritten as 

-2W (Si -Sj+ rn- 7,1-1), (6) 

where Si and rji are the spin and the pseudo spin oper- 
ators, respectively. The components of the pseudo spin 
operator are defined by 

V? = ( — 1 ) * c| T c+ x , rj7 = (-l) l c 8i c, T , 7/f = -(n iT +n i4 .-l) 

(7) 

Now, we introduce bond operators corresponding to 
the bonding and the anti-bonding modes, 

4v = 75(4+cU 4. = ^(4 - c U (8) 



The two electron states are described as B] irr A]- = 

cJg.cJ-0.. These operators on the same bond satisfy the 
anticommutation relations: 

{A ija ,A\ j(Tl } = {B ljai Bj ja ,} = <W, otherwise = 0. 

The density operators for the bond operators are given 
as 

n A a = Aj<T A H<y = \{ n i° + + T ijv)> ( 9 ) 

riBa = B\ ja Bij a = i(n itT + n, J(J - T ija ). (10) 

Since we restrict our attention only on the neighboring 
two sites we drop these indices from the operators 
defined above. 

As a trial state, we consider the following wave func- 
tion, 

\^ a )^A\ 2a Al 3s ---A[_ 1La Al M \0), (11) 

where |0) denotes a vacuum and L is the number of sites. 
This state is regarded as a Neel ordering of the bond- 
located spins, so that we call this bond Neel (BN) state. 
There is two- fold degeneracy given by and |^). In 
order to construct a model with the exact ground state, 
the local Hamiltonian hij = hij a should be decom- 
posed by the projection operators 1 — n AtT and riBa (see 
table I) in the following form, 

hij — £o = -Va^I - 71At)(1 - Ha|) + AssnsiriBj. 

+ a ab{(! - n A -[)n B i + n BT (l - n Al )}, (12) 

where eo is the ground-state energy par bond. According 
to the argument given in sec. 2, for the BN ground state, 
the parameters should be chosen as 

^AAi ^AB I ^BB > 0. (13) 

Comparing eqs. (5) and (12), the relations among the 
parameters are obtained as 

VI = |, Vj|=W, X = t — W. (14) 

The parameters in eq. (12) are identified as follows, 

\ AA = ^-W + 2t 7 (15) 

*AB = -f + W, (16) 

Abb = |+3V^-2i. (17) 

From eqs. (13), (15), (16) and (17), we obtain the pa- 
rameter space of the BN ground state as shown in fig. 2. 
Note that the BN state appears only for t > region. 





|0} Al\0) Bl\0) BlAl\0) 




10 1 
1 1 



Table I. Projection operators to construct the Hamiltonian on 
the one-dimensional chain with the bond Neel (BN) ground state. 
We can also discuss the parameter region of the ferromagnetic 
(FM) and phase-separated (PS) states. 
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A BB = U/2t 
Fig. 2. Phase diagrams of the generalized Hubbard chain (5) in 
the U/2t-W/t parameter space with t > 0. 10 ' 15 The other pa- 
rameters are set as X = t - W, Vj| = W and Vj_ = U/2. The 
shaded regions labelled by BN, FM and PS denote bond-Neel, 
ferromagnetic and phase-separated states, respectively. 



The property of the BN state can be investigated based 
on the matrix-product method. According to rcf. 10, 
both charge-charge and spin-spin correlation functions 
vanish except for those of the nearest sites which in- 
dicates the existence of the charge and the spin gaps. 
On the other hand, the bond-located spin correlation 
exhibits a long range order. We can also calculate ele- 
mentally excitation spectrum using the matrix-product 
method as a variational approach. 15 

In the present one-dimensional model at half-filling, 
we can discuss not only the BN state but also the fer- 
romagnetic (FM) and the phase-separated (PS) states. 
As shown in table I, the last term of eq. (12) stabilizes 
the fully polarized FM state for X^ B < 0. Similarly, the 
PS state where the system is separated into a domain of 
doubly occupied sites and a vacuum, is stabilized when 
^AA + Abb < 0, neglecting the surface energy. As shown 
in fig. 2, the FM and the PS states appear in the U/2t- 
W/t parameter space symmetrically in the positive- and 
in the negative-C/ regions, respectively. This is consis- 
tent with the fact that the W term is the ferromagnetic 
exchange interactions of the spins and the pseudo spins 
(6), and the PS state is regarded as the FM state of the 
pseudo-spin space. The condition W/t > 1/2 for the FM 
and the PS phases is not clearly obtained in the present 
argument. This will be discussed in Sec. 5. 

The phase boundary of the BN and the FM states 
^Ab = corresponds to the SU(2) symmetry in the spins 
Vj| = Vj_, so that the ground state is highly degenerate. 
The system undergoes a first-order phase transition at 
this level-crossing point. When W/t = 1 (X = 0), there 
is the particle-hole symmetry. At (U/2t,W/t) = (— 1, 1), 
the system has the SU(2) symmetry in the pseudo-spin 
space, so that the BN, the PS and the 77-paring states 
are degenerate. The other lines which separate shaded 
and non-shaded regions in fig. 2 do not necessarily mean 
phase boundaries. 

4. Kagome lattice 

We consider the generalized Hubbard model including 
three site terms on the Kagome lattice at 1/3-filling with 
zero-magnetic field, H = J2(ijk)a hijka, where the sum- 



Masaaki NAKAMURA and Kazuhito Itoh 3 

mation (ijk) is taken in each unit trimer as shown in 
% 1, 

+ W '{Tij a Tjka + TjkoTkia + TkiaTija) 

+X'(Ti ja n ks + T jka n ti y + T k i a nja), (18) 

where hij a is the local bond Hamiltonian defined in 
eq. (5) for z = 4. The three site interactions (W, X') 
have the similar structure to those of the t-J model. 

We define the following one-electron plaquette opera- 
tors corresponding to a zero momentum mode, left and 
right chirality modes, 

Ajka = ^(4 + 4+ c L)> (is) 

Bhk* = 73(4+^4+^U (20) 

c\ iha - ^(4+" 2 4+ wc U (21) 

where ui = e l27r / 3 . These operators on the same trimer 
satisfy the anticommutation relations: 

{Aijka, ^ijfccr'} = {Bi ]ka , Bjj ka ,} = {Cijka,Cj ]kcrl } = 5 a(7 i 

and otherwise = 0. Note that A\ jk(T \0), B\ jka \Q) and 

C\- ka |0) are chosen as eigenstates of density, hopping and 
current operators: 

Nijka = n ia +n ja +n ka , (22) 

Tijka — ^ija + ^jka ~t~ ^kia : (^3) 
Jijkcr — Jija ~1~ Jjka ~t~ Jkicr 1 (^4) 

where Jij a = i{cy ia Cj a — H.c). The density operators for 
the plaquette operators are 

n A <j = ^{N ijka +T ijk(r ), (25) 

«Bcr = -(2Nij ka — Tij ka — V3Jijkcr), (26) 
nca = g(2A r lJ fc<r - Tijka + VSJijka)- (27) 

Now we consider the Neel ordered state given by AT . fe(T 
on the honeycomb lattice as shown in fig. 1. Then the 
local Hamiltonian hij k = hij kcr corresponding to this 
state consists of projection operators 1 — UAa, »Bo and 
nca- In order to eliminate the current term from the 
Hamiltonian, we construct hij k in the following way, 

hij k - £0 = ^AAi 1 - n A|)(l - n A[) 

+Aab {(1 - n A] ){n Bi + n cl ) + (n BT + n CT )(l - n 
+Abb(?*bt + n C]){nB[ + nci), 

where the parameters should satisfy the condition of 
eq. (13). Comparing cqs. (18) and (28), the relations 
among the parameters are obtained as 

V x = —, V\\=W = W\ X = X' = t-2W 1 
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Fig. 3. Phase diagram of the generalized Hubbard model on the 
Kagome lattice (18), in the U /2t-W/t parameter space with Vj| = 
W = W, V ± = U/2, X = X' = t - 2W and t > 0. The shaded 
region labelled by PN denotes the plaquette Neel state. 



and parameters in cq. (28) are identified as 



A 



AA 



^AB 



U 

--4W 



■4t, 



= -—+4W-t, 



A 



BB 



= —+5W-2t. 



(29) 
(30) 
(31) 



Thus the parameter region of the exact ground state with 
the plaquette long-range order is obtained as shown in 
fig. 3. 

We can extend the parameter space of the same exact 
ground state, introducing more parameters in eq. (28), 
such as A BB (ri BT + n C ]){n B i + nc{) — ► ^bb { n B] n B[ + 
n CT n ci) + ^ss( ri --BT n Ci + n CT n -Bl)- Moreover, it is also 
possible to construct models for chiral states given by 
B\- k(J and Cjj ka , for the ferromagnetic state at half- 
filling, and for two-electron plaquette states at 2/3- filling 
such as 



ijka ijka 



t J 



(c T 

V 1(7 



L ja L ko 



(32) 



considered do not include hopping terms between differ- 
ent spin sectors. 

Application of this method to other lattices can be 
done straightforwardly. The bases of plaquette operators 
should be chosen reflecting symmetry of lattices. For ex- 
ample, the unit plaquettes of the Checkerboard and the 
Pyrochlore lattices consist of four sites, but the bases of 
the plaquette operators are chosen in different way. Simi- 
larly, we can also construct Hciscnberg-typc models with 
plaqucttc-ordercd ground states in these lattice systems. 
These ground states have two-fold degeneracy, but the 
uniqueness has not been proven yet. 

In this paper, we have decomposed particular Hamil- 
tonians into the projection operators. However, we need 
more general treatment to decompose the arbitrary given 
Hamiltonians. In the present case, one projection opera- 
tor gives positive value for more than two states as sum- 
marized in table. I. To remove this overlap, we should 
redefine the projection operators. For example, in one- 
dimensional case, we should choose three operators to 
treat arbitrary bond Hamiltonians as shown in table II. 
Then the number of free parameters becomes six. On the 
other hand, in sec. 3, we have introduced three param- 
eters. Among these three extra parameters, the hidden 
conditions for the ground states t > and W/t > 1/2 
are included. Moreover, when we add the three site term 
TijfjnisTijg to the model, these hidden conditions are 
essential to determine the regions of the BN ground 
state. The detail of this argument will be published else- 
where. 15 
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5. Summary and discussion 

We have discussed a method to construct general- 
ized Hubbard Hamiltonians with exact plaquettc-ordcrcd 
ground states in arbitrary dimensions. The correspond- 
ing lattices have bipartite structure in terms of corner 
sharing unit plaquettes. We have applied this method to 
the one-dimensional chain and the Kagome lattice, and 
obtained parameter regions of the exact ground states. 

In one dimension, the BN state corresponds to the 
staggered dimer states in the spin- 1/2 two-leg ladder 
model with four spin exchanges. 16 Comparing with the 
way to decompose the Hamiltonian to the projection op- 
erators in ref. 16, the present argument is quite simple. 
This is because the generalized Hubbard models we have 





|0) Al\0) Bl\0) BlAl\0) 


(1 - nA CT )(l - n Ba ) 
(1 - n Aa )n Bcr 

nAtrlT-Bcr 


10 
1 
1 



Table II. A different set of projection operators for the one- 
dimensional chain with the BN ground state. 



1) J. Hubbard: Proc. R. Soc. London A 276 (1963) 238; J. 
Kanamori: Prog. Theor. Phys. 30 (1963) 275; M. C. Gutzwiller: 
Phys. Rev. Lett. 10 (1963) 159. 

2) D. K. Campbell, J. T. Gammel and E. W. Loh, Jr.: Phys. Rev. 
B 38 (1988) 12043; 42 (1990) 475. 

3) R. Strack and D. Vollhardt: Phys. Rev. Lett. 70 (1993) 2637. 

4) R. Strack and D. Vollhardt: Phys. Rev. Lett. 72 (1994) 3425. 

5) L. Arrachea and A. A. Aligia: Phys. Rev. Lett. 73 (1994) 2240. 

6) J. de Boer, V.E. Korepin and A. Schadschneider: Phys. Rev. 
Lett. 74 (1995) 789. 

7) J. de Boer and A. Schadschneider: Phys. Rev. Lett. 75 (1995) 
4298. 

8) A. Montorsi and D. K. Campbell: Phys. Rev. B 53 (1996) 5153 

9) M. Kollar, R. Strack and D. Vollhardt: Phys. Rev. B 53 (1996) 
9225. 

10) K. Itoh, M. Nakamura and N. Muramoto: J. Phys. Soc. Jpn. 
70 (2001) 1202. 

11) K. Itoh: J. Phys. Soc. Jpn. 68 (1999) 322. 

12) A. Tanaka and H. Ueda: Phys. Rev. Lett. 90 (2003) 067204. 

13) Y. Imai, N. Kawakami and H. Tsunetsugu: Phys. Rev. B 68 
(2003) 195103. 

14) for example, O. A. Petrenko and D. McK. Paul: Phys. Rev. B 
63 (2001) 024409. 

15) M. Nakamura, T. Okano and K. Itoh: unpublished. 

16) A. K. Kolezhuk and H.-J. Mikcska: Phys. Rev. Lett. 80 (1998) 
2709; Int. J. Mod. Phys. B 12 (1998) 2325. 



